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O ■ Abstract 

We study the properties of time evolution of the — system 
in spectral formulation. Within the one-pole model we find the exact 
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form of the diagonal matrix elements of the effective Hamiltonian for 
(~| \ this system. It appears that, contrary to the Lee-Oehme-Yang (LOY) 

result, these exact diagonal matrix elements are different if the total 
system is CPT-invariant but CP-noninvariant. 



1 Introduction 

This paper has been inspired by the results presented in [l] and [2]. Pa- 
per 111 analyses the problem of equality of particle and antiparticle masses, 
whereas |2| describes an exactly solvable model of the particle-antiparticle 
system - in this particular case — K^. The most important properties 
of antiparticles follow from the CPT symmetry. This symmetry, also known 
as the CPT theorem p], determines the properties of the transition ampli- 
tudes under the action of the product of C, P and T transformations (charge 
conjugation, space inversion and time reversal, respectively). According to 
the CPT theorem, the transition amplitudes describing any physical process 
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must be CPT-invariant. From this we conclude that the full Hamiltonian H 
of the system under consideration must be invariant under the product of 
the C, V and T operators. Another conclusion that can be drawn here is 
that stable particles and their antiparticles must have the same mass. This 
property of the particle-antiparticle pair is true for stable particles and the 
same is usually assumed of unstable particles (e.g. and K^). Such an 
extension of a law true for stable particles to unstable particles is questioned 
in [l]. The reason for the widespread belief that this extension is correct is 
most probably the Lee, Oehme and Yang (LOY) approximation and the con- 
clusions which follow from it — and more specifically, the properties of the 
effective Hamiltonian, Hloy, governing the time evolution in the subspace 
7i||. In our case TC\\ is the subspace of the total Hilbert space of states H, 
spanned by state vectors of K^, mesons. 

Following the LOY approach, a nonhermitian Hamiltonian ify is usually 
used to study the properties of the particle-antiparticle unstable system |H- 

Hii^M- (1) 

where 

M = M+ , F = r+ (2) 

are (2 x 2) matrices acting in 7i\\. The M-matrix is called the mass matrix 
and F is the decay matrix. Lee, Oehme and Yang derived their approx- 
imate effective Hamiltonian H\\ = Hloy by adapting the one-dimensional 
Weisskopf-Wigner (WW) method to the two-dimensional case correspond- 
ing to the neutral kaon system. Almost all properties of this system can be 
described by solving the Schrodinger-like equation 

t^^\t,t)\\=H\^\ij;t)\i, (t>to>-oo) (3) 
(where we have used h = c = 1). The initial conditions for Eq. ^ are 

II |^;t = to)|| 11=1, |V';to = 0)|| =0, (4) 

where lip] t = to)\\ belongs to H\\ {H\\ C H) and Tiy is spanned by orthonormal 
neutral kaons states: |ir°) = |1), |i?°) = |2). Thus n\\ = PH, where 

P=|l)(l| + |2)(2|. (5) 
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According to the standard result of the LOY approach, in a CPT invariant 
system, i.e. when 

ene-^ = h, (e) 

(where 6 = CPT) we have 

.LOY _ uLOY 



h-r = (7) 



and 



= M^OY, (8) 
where: ^fj^^ = ^(^j/^^) ^{z) denotes the real part of a complex 
number z (then '^{z) is the imaginary part of z), and hjj^^ = (jl-f^Loylj) 
(J = l,2). 

The universal properties of the two particles subsystem described by the 
H fulfilling the condition (jH)), may be investigated by the use of the matrix el- 
ements of the exact evolution operator for 7i|| instead of the approximate one 
used in the LOY theory. This exact evolution operator, U\\{t), can be writ- 
ten as follows U\\{t) = PU{t)P, where U{t) = e~^^^ is the exact evolution 
operator acting in the total state space Ti. 

Assuming that the CPT symmetry is conserved in the system under con- 
siderations one finds that the matrix elements 

A^k{t) = m\\{tm^m{tm (j,^ = i,2), (q) 

of the exact U\\{t), obey 

All(t) = A22{t). (10) 

General conclusions concerning the properties of the difference of the 
diagonal matrix elements {hu — /i22) of the exact H\\, (which can in general 
depend on time t [E]), where 

h,k = m\\\^ {3,k = 1,2), (11) 

can be drawn by analyzing the following expression derived in pj for CPT 
invariant systems 

In [l] it is shown that 

hu{t) - h22it) ^ 0, (13) 
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when (in)) holds and 

[CP,F]^0, (14) 

that is in the exact quantum theory the difference (/iii(t) — /i22(^)) cannot 
be equal to zero with CPT conserved and CP violated. In Section 3 we 
will consider this relation in the context of an exactly solvable model. This 
problem is important because realistic calculations are carried out with the 
use of simplified and approximate models. Not all of them conform to the 
requirements of the exact (not approximate) quantum theory. 

The aim of this paper is to calculate the difference of the diagonal ma- 
trix elements of the effective Hamiltonian, (fT2|l . in a CPT invariant and CP 
noninvariant system for the approximate model analyzed in P], that is in the 
case of the one - pole model based on the Breit - Wigner ansatz, i.e. the 
same model as used in Lee, Oehme and Yang theory. 

The paper is organized as follows. In Section 2 we review briefly the 
spectral formulation for the neutral kaon system and a model described in [2]: 
one pole approximation. Section 3 investigates the diagonal matrix elements 
of the effective Hamiltonian and their difference in the CPT invariant and 
CP noninvariant system in the case of the one - pole model. In Section 
4 we present our conclusions and we estimate the numerical result of the 
investigated difference for the K'^ — system. Appendix A contains the 
relevant integrals and derivatives used in Section 3. In Appendix B we give 
the exact formulae for expressions appearing in Section 3. In Appendix C 
we find coefficients which were used in Section 3. 

2 The model: one pole approximation 

While describing the two and three pion decay we are mostly interested in 
the \Ks) and \Ki) superpositions of \K^) and \K^). These states correspond 
to the physical \Ks) and \Ki) neutral kaon states [2l fTT| IT2] 

\Ks)=p\K'^) + q\K'^), \KL)=p\K^)-q\K''). (15) 

We assume that these physical states are the initial physical states of the 
CPT - invariant system, i.e. at the instant of creation of neutral kaons. We 
have 

{Ks\Ks) = {Kl\Kl) = \pf + \qf = 1 (16) 
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{Ks\Kl) = {Kl\Ks) = bp - \q\' Ak ^ 0. (17) 

The time evolution of K'^ and can be concisely presented in the fol- 
lowing way: 

where Ka = K^, and H is the full hermitian Hamiltonian and Q = I— P, 

AK^K,{t) = {K^le-'^^'lKp) = {K^\Kp{t)). (19) 

Let us notice that amplitudes AK^Kfjif) for Ka^ Kp = K^, correspond 
to the previously defined amplitudes^ Aji.{t), where j,k = 1,2, Q. Conse- 
quently we may write 



ARORoit) = 




e" 


-itH 


\K^) = 


(l|e- 


-itH\ 


|1) 


= Ai(t), 


ARORoit) = 






-itH 


\K^) = 


(l|e- 


-itH\ 


|2) 


= Al2(t), 


AR0R0{t) = 




\e' 


-itH 


\K^) = 


(2|e- 


-itH\ 


|1) 


= ^2l(t), 


AR0R0{t) = 




\e' 


-itH 


|K°) = 


{2\e 


-itH\ 


|2) 


= ^22 (t). 



(20) 

Using the spectral formalism we can write unstable states |A) as 

|A) = E^A(g)|g) (21) 

and then |A(t)) as 

|A(t))'=^-^*^|A) = El9(t))^A(g), (22) 

where \q{t)) = e~**^|g) and vectors |g) form a complete set of eigenvectors of 
the hermitian, quantum-mechanical Hamiltonian H and uj\{q) = {q\^)- If the 
continuous eigenvalue is denoted by m, we can define the survival amplitude 
A{t) (or the transition amplitude in the case of ^ ) in the following 
way: 

A{t)= J dme-''^'p{m), (23) 

Spec{H) 



^Amplitudes AKaKa{t), etc., correspond to PKOK°{t), etc. respectively used in (2] 
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where the integral extends over the whole spectrum of the Hamiltonian and 
density p(m) is defined as follows 

p(m) = \uJx{ni)\^, (24) 

where ujx{m) = (m|A). 

The above formalism may be applied to \Ks) and \Kl) by introducing 
a hermitian Hamiltonian with a continuous spectrum of decay products la- 
belled by a, P etc., 

H\(f)a{rn)) = m |</)a(m)), (0/3(m')|0a(m)) = SapSim' - m). (25) 

In accordance with formula (f22|l the unstable states Ks and Kl may now be 
written as a superposition of the eigenkets |0o(m)), 

\Ks) = rfm ^u;5,^(m)|0„(m)), (26) 

roo 

\Kl)= rfm 5]^i,^(m)|0;3(m)). (27) 

Jo a 



Thus 



\Ksif)) = e-^'^'lKs) = / dm ^ cu5,,(m)e-**^|0„(m)). (28) 

Jo 

\KL{t)) = e-**^|Ki) = / dm Y.ULAm)e-''''\Mrn))- (29) 

"'0 a 



Using (f28|) and ()29|1 we can write 

{Ks\Ks{t))= Hdm Y.\^sA^)?^~''^' = ^KsKsit), 
Jo Q 

{KL\KL{t)) = TdmY. |u;^,^(m)|V^-* ^x,x,(t), 

{Ks\KL{t)) = r dm Y.ujl^{m)ujL,,{m)e-'^'''^ AKsK^t). 
{KL\Ks{t)) = r dm Y.ulAm)uJsAT^)e-''^' = ^K,Ks{t). (30) 



From (fT^ we can obtain 

2p 

and 



K') = i:i\Ks) + \K,)) (31) 



= ^^{\Ks) - \K,)). (32) 

Now, using formulae (jSHj), (f^ - (|r?IHl . we can express AxoKo{t) etc. in 
terms of quantities describing physical states, that is through the amplitudes 
AKsKs{i)i ^KLKsit) etc. (see eg. [21) 

Akok<^ (t) = ^ [AksKs (t) - Ak,k, (t) - AksK, it) + Ak.Ks it)] , (33) 

Aroro (t) = ^ [AksKs (t) - Ak,k, (t) + AksK, {t) - Ak.Ks {t)] ■ (34) 
One can also find that 

AksKs + Ak^K^ = 2 {\p\^AKOKo{t) + \q\^AKOKo{t) ). (35) 

Assuming (jH} and using (fT6|) . (f35|) we get 

AKoMt) = Aj^oRoit) = i [AKsKsit) + AK.KAt)) ■ (36) 

It follows from ((HH), JSl and JSSl) - dSEl) that the probabilities AxoKoit) etc. 
can be written in the following way: 



Aft'O^o(t) = A^o^o(t) = / dm pft:o^o(m)e 

Jo 

= dmJ2{\usArn)\' + M\m)}e-^^' (37) 

Z Jo Q 



("OO 

A^oi^o(t) = / dm pf(OKoim)e~'"^^ 
Jo 

1 f 



-a;*^^(m)c^i,^(m)+u;2,/3("^)^s,/3(m)|e (38) 
7 



AjiOKoit) = I dm pKOKoirn)e 
1 /""^ f 



+ujlp{m)ujLA^) - ^UM^sA^) e (39) 



The Breit-Wigner ansatz ^ 

r 1 

pBw{m) = —- , = kMP (40) 

/TT (m — mo) + Y 

leads to the well known exponential decay law which follows from the survival 
amplitude 

/OO 111 
dm e-'^^'pBwim) = e-'^^^^e-^ri*!. (41) 
-00 

(Note that the existence of the lower bound for the energy (mass) induces 
non-exponential corrections to the decay law and to the survival amplitude 
PT|l — see |2| ). It is reasonable to assume a suitable form for u!s,/3 and u!l,i3- 
More specifically, we use [2J 

^s,/3 {m) = \— — , (42) 

V Ztc m — ms + t-^ 

^lA^) = \ TT- —TT (43) 

V ZTT m — m^ + ^ 2 

where As^p and Ai^p are the decay (transition) amplitudes. It is convenient 
to use the following definitions: 

75 = 7l = Am = mL-ms (44) 

s^Y.\Ma\^ l^J2\^lA' (45) 

a a 

R^Y. mh^L,.), I = E ^>(^5,.^L,.). (46) 
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In the one-pole approximation ((l^ . A^Of^olt) can be conveniently writ- 
ten as 



Collecting only exponential terms in one obtains an expression analogous 
to the WW approximation 



Here NxoKoit) denotes all non-oscillatory terms present in the integrals (|T7jl . 

3 Diagonal matrix elements of the effective 
Hamiltonian 

This section constitutes the main part of the paper. Using the decomposition 
of type (jiSll and the one-pole ansatz (j42| . (|43ll . we find the difference (fT3| . 
which is now formulated for the — system. It can be written as follows 




(47) 



A^OKoit) = A^QKoit) = - (e 



hllit) - h22{t) 



(49) 



Y{t) 



where 




(50) 



and 



Y{t) = AKOKoit)AjioKo{t) - Aj^oKoit)AjioKoit). 

To calculate p8|) . (f39|) we use the following relations [2] 



(51) 




(52) 
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and 



roo , 

JO p 



m 



~imt 



dm 



aiw? + him + ci 



-imt 



[{m - msY + 7|] [(m - mLY + Tl] 



TT 



75 



75 



+75C, (^-J«(75i,-^)) +K^^\ist)^ 



IL ' 



(53) 



where ai,6i,Ci and Cj,D'j,Fj are defined in Appendix B and J'^^\a^'q)^ 
J^^\a,r]), K^'^^a), K^^\a) in Appendix A. 

Using the above mentioned formulae from Appendixes A and B (without 
any additional simplifications and approximations) we get, for example 



it) = 



TT 



S7rp*q 



_^-im.Lt^-jLt 



l-kr 



1 + ks 

+NKOKo{t), 



where 



75 V 



-2i^sCi + D'j-F'j 



(54) 



(55) 



(56) 



and Nj^oKO (t) is the non-oscillatory term containing the exponential integral 
function Ei and it has following form 

Smp*q [ 
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-2iCi + —{-D'j + F[) 

Is 



2iCi + —{D'j-F[) 



+ ^/lSlL 



(57) 



Using the expression (jHT|l for the derivative of Ei (Appendix A) we can 
find the derivatives which will be necessary for the following calculations, for 
example 



M^o^o(t) 

dt 



1 + TT 



-imst --^St I 



-inis - 75(1 + ks) 



^^-im,t^-^,t 1^^^^ _ ^^^^ ^ 
+AiV^Oi?o(t), 

where ANj^oxoit) is defined as follows 



(58) 



AiV^o^o(t) = -l—le-'"'^'e-''^'E,{^st + tmst)(-zms-^s{l + ks)] 
+e-'^^'e-^^'E,{^Lt + iniLt) [imL - 7l(1 + k^)^ 
+e-'^''e^''E,{-^st + inist) {^nis - 7s + 
+v/7i7l(-2^7sC/ - + F;)^ 
+e-*'"^*e^^*Ei(-7Lt + iniLt) (^irriL - 7l + 
+V7i7Z(2mC/ + /^;-F;)^|. 



(59) 



There are similar expressions for Ap^oj^o{t), Nj^oKo{t), ^ii^MlW.^ ANj^oKo{t). 
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The states \Ki) and \Ks) are superpositions of \K^) and \K^) (JSIJ, 
(jn2I))- The lifetimes of the \Kl) and \Ks) particles may be denoted by tl 
and r^, respectively, = = (5.17 ±0.04) ■ 10~^s being much longer than 
Ts = ^ = (0.8935 ± 0.0008) ■ lO-^^s 

Below we calculate the difference (|i9|) for t tl 



hii{t ~ Tl) - h22{t ~ Tl) 



X{t ~ Tl) 
Y{t ~ tl) 



(60) 



When we consider only the long living states \Kl) then we may drop all 
the terms containing e~'^^^\t^T-L they are negligible in comparison with the 
elements involving the factor e~'^^*|j^T-i- We also drop all the non-oscillatory 
terms Nxoxoit), iV^o;^o(t), iV;^o^o(t) dSIl) present in A^OKoit) ((HH)), v4j^o^o(t) 
and A;^o^o(t) (fMj) . because they are extremally small in the region of time 
t ^ Tl [21 • Similarly, because of the properties of the exponential integral 



function Ei, we can drop terms like ANj^op^o in 



dA 



at 



and AN 



RORO 



jnni) 



m 



at 



(f58|) . This conclusion follows from the asymptotic properties of the 
Ei function (|80| and the fact that ANj^oko, ANj^oxo only contain expressions 
proportional to E^. 

We may now calculate the products AK<>Ko{t)-Aj^oKo{t), y4^o^o(t)-A^oxo(t), 



■(if:) ■ Aj^oKo{t), that after the use of the above 



mentioned properties of Nxoxoit), AA'"^o^o(t) and performing some algebraic 
transformations, lead to the following form of the difference (|60|l : 



huit ~ Tl) - h22(t ~ Tl)) 



7r2 + 27r + l/ V 



7^0, 



(61) 



where 

z = Mp\"\q 



2| |2 7r2 + 277 + 1 
47]-2 



+AiCj{D'j - f;) 



7^0 



(62) 



W = 2{ -CiniL + Dj-Fj]+i 



-4CjjL + —{-D'j + F[] 



+ 0. (63) 
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4 Final remarks 

Our results lead to the conclusion that in a CPT invariant and CP noninvari- 
ant system in the case of the exactly solvable one-pole model, the diagonal 
matrix elements do not have to be equal. In the general case the diagonal 
elements depend on time and their difference, for example at t ~ r^, is differ- 
ent from zero. This has been clearly demonstrated in the last Section: Z and 
W in (|6T| are different from zero, so the difference (/iii(t) — /i22(^))|t~Ti, 7^ 0. 
From the above observation a conclusion of major importance can be drawn, 
namely that the measurement of the difference (/iii(t) — h22{t)) should not 
be used for designing CPT invariance tests. This runs counter to the general 
conclusions following from the Lee, Oehme and Yang theory. 

A detailed analysis of hjk{t), {j,k = 1,2) shows that the non-oscillatory 
elements Na^pit) , ANa^^lt) (where a, = K^,T{^) is the source of the non- 
zero difference (/iii(t) — /i22(^)) in the model considered. It is not difficult to 
verify that dropping all the terms of Na^/sit) , ANa^^it) type in the formula 
for {hn{t) - h22{t)) gives (/i?r(t) - = 0, where h°f{t), (j = 1,2), 

stands for hjj{t) without the non-oscillatory terms. 

To obtain the numerical estimate the real and imaginary parts of 
{hiiit ~ tl) — h22(t ~ tl)) it is necessary to put experimentally obtained val- 
ues of ■ms,mL,'ys,'jL, etc., into (|6T|) - (|63|) . According to the literature OEI, 
if the total system is CPT-invariant but CP-noninvariant then we have (see, 
eg., M) 

1 + e 1-e 

and hence we get 

AK = 2-^{e) (65) 

where |£:| ^ 10~^ [HIIZI- Putting experimental values jTH] 

Am = (3.489 ± 0.008) x 10~^^MeV, (66) 
m5 ~ mi ~ maverage = (497.648 ± 0.022)MeV" (67) 

and Tl, ts, h into expressions (jHT|l - for the neutral kaon system we can 
obtain the following estimations 

3? (/iii(t ~ Tl) - h22{t ~ Tl)) - -4.771 x IQ-^^MeV (68) 

13 



(64) 



and 



$5 {hn{t ~ tl) - h22it ~ tl)) ~ 7.283 x IQ-'^MeV. (69) 

So, the difference of the diagonal matrix elements of the effective Hamiltonian 
for the — system in one pole approximation is different from zero. 
According to our evaluation 



|3?(/iii(t ~ Tl) - h22{t ~ tl))\ _ \mKO - mKo\ 



10-^\ (70) 



^^average "^average 



Recent experiments give ^"^-^^ '"-^'-'"^ < 10 Il5j. So our estimation (ITOl) does 

'^average. ^ ^ ■ ^ 

not contradict the experimental results. 

Deviations from the LOY result estimated in [2] have the order of magni- 
tude ^. These estimations refer to amplitudes A^^j^o and A^^Ko- However, 
these estimations could not be directly transformed into the calculation of 
the difference (/iii(r) — /i22(T)), because the difference depends not only on 
amplitudes of type Aj^^j^^, but also on their derivatives (see: relations (|i9|) - 

(jnij). There are products of type ^^^^^^^Aj^oxoit) in the numerator of the 
expression whereas there aren't any derivatives in the denominator of 

this expression. What is more, there is the difference of expressions of type 
^"^^af ° Akoro [t) in the numerator of (j49| . So, it can hardly be expected, 
that the order of deviations from the LOY result of the relatively complicated 
expression will be the same as the order of corrections to the LOY result 
of one of the following expressions: ^XoXo ^^'^ ^/?o^o- 

If estimation (|70|) is compared with a similar one obtained in [l], [H] - 
[TS] one can see that the numerical value of our estimation is much larger 
than the value of mentioned estimations. It is because the estimations given 
in the mentioned papers were obtained using a different method for the Lee- 
Fridrichs model [T9] . 

The results {hii{t) - /i22(^)) 7^ and ((Ml), (EHl) and (HOI) seem to be very 
important as they have been obtained within the exactly solvable one-pole 
model based on the Breit-Wigner ansatz, i.e. the same model as used by Lee, 
Oehme and Yang. 

As the final remark it should also be noted that the real parts of the 
diagonal matrix elements of the mass matrix H\\, hu and h22, are considered 
in the literature as masses of unstable particles |1), |2) (e.g., mesons Kq, Kq). 
The interpretation of the diagonal matrix elements of H\\ (t = 0) = PHP is 



14 



obvious (see P^). They have the dimension of the energy (that is, the mass) 
and hjj{0) = {j\H\j), (j = 1,2). So their interpretation as masses of particle 
"1" and its antiparticle "2" at the instant t=0 seems to be justified. Note 
that H\\ has the following form ([IHl, [HI - [HI) 



The diagonal matrix elements of the operator V||(t), i.e. Vjj{t) = {j\V\\{t)\j) , 
also have the dimension of the energy and in general they depend on time. 
So, the problem seems to be open: we can treat the matrix elements of the 
operator V||(t) as a time-dependent correction to the energy or mass. 
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A Appendix 

This appendix contains the relevant integrals and derivatives used in Section 
3. 

Integrals K^"^ (a) and J*^"^ (a, i]) are defined as follows [2 [201 



PHP+V\\{t) 



(71) 



that is 



H\\{Q) + V\\{t). 



(72) 




(73) 




(74) 



If we assume a = {"is/ hi) and rj = ( 



IS/L 



), for n = we get 




15 



-IS/L 



(75) 



IS/L 



'^S/L 



dy 



2/2 + 1 



-iis/L^y 



1 



.^sgn{-'^)e-'^s/Lt 



IS/L 



+e^s/^*E,(-7s/Lt) 



-is/ht 



IS/L , 



ISIL 



7S/L ' 



(76) 



stands for the 



where Ei is the exponential integral function and sgn i 
sign of (-^). 

° ^ 7S/L ^ 

Any other integral of K^'^\a) or J^"\a, r]) for n > can be obtained from 
(f73|l or (f7i|l by differentiating (f73|) or (fTHl with respect to a and using the 
Fourier identity in (f7i|l 0. For n = 1 we have 



dy- 



X 



+ 
2 



+e^s/^*E,(-7s/Lt) 



(77) 



75/L 



X 



2/2 + 1 



-i^s/Lty 



1 

"2 



75/L 
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Ei Js/Lt[l - 



+e^s/^*E,(-7s/Lt) 



IS/L 

1s/l' 



(78) 



The exponential integral function Ei is defined in the following way [2], [201 



E,i±xy) = ±e^^y I 
Jo 



dt , 

vTt 



> 0, a; > 0. 



(79) 



We can use the very convenient asymptotic properties of Ei given in [21] 

E,(0) = -oo, 
Ei{oo) = oo, 
E,(-oo) = 0, 
Eiiioo) = in, 

E,{-ioo) = -i-K. (80) 

These properties of Ei have been used to obtain the final result (|IITll - 

In our calculations we have also used the formula for the derivative of Ei. 
Its final, general form is given below 

dEi{±xy) 1 



dx 



— e 

X 



±xy 



(81) 



B Appendix 



In this Appendix we collect from [2] the coefficients ai,6i,ci and Ci,D'j,Fj 
which were used in Section 3. 

The calculations will be clearer if we write the sum of the product 
J2i3^^s,i3^L,i3 in the same way in which spectral functions defined by (|i2|l . 
were used earlier in |E1 



BW 



7r[(m - msY + 75] [("^ - + ll] 



X 



x{(aj?m^ + + cr) + i^ajm^ + bjm + c/)} 



(82) 
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where 

aj = J, bi = (75 - 7L)i? - {ms + mi)! 
ci = {iLfns - 1smL)R + {rnsmi + 757^)/ 

similar formulae may be found for a^j, 6r, c/j, where 

R = ,r^^{lsS + lLL), 



Am 

-(75^ -7l^), 



and 



Tcms 

L = l + ^ + 0((7i/m^)2). 



Equations ()86|1 result from performing the following integration 

/•OO /"OO 

/ dm y2\^s,a\'^ = dm y2\^L,i3\'^ = 1- 

This integral follows from the initial conditions defined by (fT5|) - (fT7jl . 
This expression is now factored 

a/m^ + bjm + c/ 



[(m - msy + 7|] [{m - mi^f + 7I] 
Cim + Di _^ Ejm + Fj 



{m — ms)"^ + 'Js {m — miY + 

which leads, as usual, to a linear system of equations which allows 
calculate the coefficients Cj,Dj,Ej, Fj 

Ej = —Ci, 

CiAm + D'j + F'j = ai 

Ciiiml + 7I) - (m| + 7I)] - 2D'j{mL - 2F[ms = h 
D'Aml + 7I) + F[{ml + 7I) + C/[mi(m| + 7I) 
-m5(m^ + 7^)] = c/ 
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In the last formula we have introduced new definitions 



Ej = —Ci, 

D'j = Dj + CjmL, F^ = Fj-CjmL. (90) 
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